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The dynamics of symmetry breaking responsible for lifting the degeneracy of the Landau levels 
(LLs) in the integer quantum Hall (QH) effect in graphene is studied in a low-energy model with 
the Coulomb interaction. The gap equation for Dirac quasiparticles is analyzed for both the lowest 
, and higher LLs, taking into account the LLs mixing. It is shown that the characteristic feature of 

the long-range Coulomb interaction is the dependence of the gap parameters on the LL index n 
' ("running" gaps). The renormalization (running) of the Fermi velocity as a function of n is also 

studied. The solutions of the gap equation reproduce correctly the experimentally observed integer 
QH plateaus in graphene in strong magnetic fields. 

PACS numbers: 73.22.Pr, 71.70.Di, 71.70.-d 



I. INTRODUCTION 



As it is well known, the low-energy dynamics of electrons in graphene [l[ is described by the Dirac equation in (2+1) 
dimensions [2j . Perhaps the most direct confirmation of the pseudorclativistic character of electron motion in graphene 
is given by the experimental observation 0, Q of the anomalous quantum Hall (QH) effect theoretically predicted in 
Refs. The anomalous QH plateaus are observed at the filling factors v = ±4(n + 1/2), where n = 0, 1, 2, . . . 

is the Landau level (LL) index. The factor 4 in the filling factor is due to a fourfold (spin and sublattice-valley) 
degeneracy of each QH state in graphene. The presence of the anomalous (from the viewpoint of more standard 
condensed matter systems) term 1/2 in the filling factor unmistakably reveals the relativistic-like character of electron 
q ■ motion in graphene \$L Icl- 111. 

The later experiments 12, [H| m strong magnetic fields (B > 20 T) observed new QH plateaus, with integer filling 
factors v = 0,±1, and ±4. The more recent experiments [3, EH discovered additional plateaus, with v = ±3 and 
v ± 1/3. While the latter corresponds to the fractional QH effect, the plateaus with v = 0,±1,±4, and v = ±3 are 
, intimately connected with a breakdown of the £/(4) symmetry of the low-energy effective quasiparticle Hamiltonian in 
graphene (connected with the spin and sublattice-valley degeneracy mentioned above) [Toj . Strictly speaking, because 
C*") , of the Zeeman effect, the f/(4) is reduced to a Uf(2) x U±(2), with U(2) s being the sublattice-valley symmetry at 
a fixed spin (s —'f or s =i). However, taking into account that the Zeeman interaction is rather weak for realistic 
l/"") \ magnetic fields, the U(A) is a good approximate symmetry guaranteeing that at weak magnetic fields only the QH 
• plateaus with the filling factors v = ±4(n + 1/2) appear. The observed new QH plateaus v = 0, ±1, ±4, and v = ±3 
occur clearly due to the electron-electron interaction leading to (quasi-)spontaneous U(A) symmetry breaking that 
removes the degeneracy of the n = and n = 1 LLs. 

For the description of the new QH p lateaus, the following two theoretical scenarios were suggested. One of them is 
the QH fcrromagnetism (QHF) [16l4l9l |. whose order parameters are the spin and valley charge densities (the dynamics 
of a Zeeman spin splitting enhancement considered in Ref. [2(| is intimately connected with the QHF). This scenario 
is related to the theory of exchange-driven spin splitting of LLs [2l| . The second one is the magnetic catalysis (MC) 
scenario, whose order parameters are excitonic condensates, responsible for the generation of the Dirac masses of 
charge carriers [22Tj25j . The essence of the magnetic catalysis, which is connected with the effective dimensional 
reduction in the dynamics of charged fermions in an external magnetic field, was revealed in Ref. [26| . It was first 
applied for a single layer of graphite in Refs. [9l. [Toj. 
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One may think that the QHF and MC order parameters should compete with each other. However, the analysis 
in an effective model with the local four-fermion interaction performed in Ref. [27j showed that these two sets of the 
order parameters necessarily coexist (this feature has been recently discussed also in Ref. [13 )• This fact strongly 
indicates that these two sets of the order parameters have a common dynamical origin, i.e., they are two sides of the 
same coin. Their simultaneous consideration qualitatively reproduces all the QH plateaus observed experimentally in 
strong magnetic fields. 

Certainly, it would be important to extend the analysis of Ref. [27| to the case of realistic long-range Coulomb 
interactions. In the present paper, the dynamics of C/(4) symmetry breaking responsible for the appearance of the QH 
plateaus with v = 0, ±1 in the lowest LL and with ±3, ±4, and ±5 in the n = ±1 LL is studied in a low-energy model 
with the Coulomb interaction. While the symmetric structure of the solutions descibing these plateaus is similar to 
that of the solutions in the model with the local four-fermion interaction [27j , there are essential qualitative differences 
between them. In particular, the long-range Coulomb interaction leads to the decrease of the gap parameters with 
increasing the LL index n ( "running" gaps) . 

Recently, the dynamics with the Coulomb interaction in the v = quantum Hall state was studied in Refs. [28l . [29| 
by utilizing different approaches than the present one. In these papers, the important role of the Landau levels mixing 
effects was also revealed. Note that these mixing effects are important both in models with long-range interactions 
and those with short-range ones, as the model considered in Ref. p7j |. 

The rest of the paper is organized as follows. In Sec. HH the general features of the model, in particular the 
structure of the order parameters, are described. The Dirac quasiparticle propagator and energy dispersion relations 
are considered in Sec. IIII1 The Schwinger-Dyson (gap) equation for the quasiparticle propagator is derived in Sec. IIVI 
In Sec. |Vl we present our numerical results, which include (i) the renormalization of the LL-dcpendent Fermi velocity 
parameter, (ii) the solutions of the gap equation at the lowest Landau level (LLL) and (iii) the solutions of the gap 
equation at the n = 1 LL, respectively The general discussion of the main results and a summary are given in Scc. lVIl 
Several appendices at the end of the paper contain some of our derivations and technical details used in the main 
text. 

Note that in most parts of the paper, we use the units with K = 1 and ks = 1. 



The low-energy quasiparticle excitations in graphene are conveniently described in terms of a four-component Dirac 
spinor \& s = {ipKAs,'4 ; KBs,'>pK'Bs,'4 , K l As) which combines the Bloch states with spin indices s =f or s =1 on the 
two different sublattices (A, B) of the hexagonal graphene lattice and with momenta near the two inequivalent points 
(K, K') at the opposite corners of the two-dimensional Brillouin zone. The free quasiparticle Hamiltonian has a 
pseudorelativistic form with the Fermi velocity vp w 10 6 m/s playing the role of the speed of light 



where r = (x,y), \I/ S = is the Dirac conjugated spinor and summation over spin is understood. In Eq. (J]), 

7 M with /i = 0,1,2 are 4x4 gamma matrices belonging to a reducible representation 7 M = T3 ® (T3, 1T2, — iri) 
where the Pauli matrices f,r act in the subspaces of the valley (K, K') and sublattice {A,B) indices, respectively. 
(For the Dirac 7-matrices, we use the same representation as in Ref. |22j, Appendix C.) The 7-matrices satisfy the 
usual anticommutation relations {7^,7"} = 2g fJ - l/ , = diag (1, — 1, — 1) , fi, v = 0,1,2. The canonical momentum 
tv = (TT x ,Tr y ) = —ifiV + eA/c includes the vector potential A corresponding to a magnetic field B±, which is the 
component of the external magnetic field B orthogonal to the graphene xy-plane. 
The Coulomb interaction term has the form 



where Uc(r) is the Coulomb potential in a magnetic field. The corresponding potential with the polarization effects 
taken into account was represented, for example, in Eq. (46) of Ref. PI|. The Hamiltonian H = Hq + He possesses 
the t/(4) symmetry discussed above. The electron chemical potential fi is introduced through adding the term —fi^ty 
in H (this term preserves the E/(4) symmetry). The Zeeman term ^jMVs^, where B = \B\ and hb = eh/(2mc) is 
the Bohr magneton, and 173 is the third Pauli matrix, breaks the t/(4) symmetry down to the U^(2) x U±(2) symmetry. 
The generators of the U^(2) x U±(2) subgroup are given by 



II. MODEL 




(1) 




(2) 
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where 7 s = i7°7 1 7 2 7 3 and P s = (1 + sa^)/2 are the projectors on spin-^ (s = +1) and spin- 4. (s = — 1) states. Note 
that U s (2) ~ U s (l) x SU S (2), where the first matrix in Eq. ^ is a generator of the U s (l) and last three matrices are 
the generators of the SU S (2). 

The order parameters that describe the breakdown of the U(A) symmetry are the same as in Ref. [27J ■ The charge 
densities (and corresponding chemical potentials), which span the QHF order parameters, are given by: 

M3 : *W= E (^L^-^U^ai), (4) 

n=K,K' a=A,B 

fl s : *V7 5 -Pa* = ^KAs^KAs + i>K Bs ^KBs ~ ^K'As^K'As ~ Bs^K' Bs , (5) 

where 

TV-(J_°r) (6) 

is the diagonal valley matrix related to the SU S (2) symmetry. The chemical potentials (13 = (jj^ — fi\)/2 and fl s are 
related to the spin and valley densities, respectively. The latter describes a charge density imbalance between the two 
valleys in the Brillouin zone (or the anomalous magnetic moment in the language of relativistic field theory) . 

Their MC cousins are connected with charge density wave (CDW) and valley polarized CDW. The corresponding 
masses are the Dirac mass and the Haldane mass [§| , 

A s : *P S * = ^ KAs lj}KAs ~ ^KBs^KBs + ^ K , As ^K'As ~ ^ K , Bs ^K'Ba , (7) 

A, : * 7 3 7 5 P S * = ^ ka Jka s - ^ KBs ^kb s - (<- As #'A s - ^ k , Bs ^k'b s ) , (8) 

respectively. 

While the generation of CDW in Eq. ([7]) breaks spontaneously the U(2) s [more precisely, its subgroup SU(2) S ], the 
generation of valley polarized CDW in Eq. (JSJ does not break this symmetry. On the other hand, while the Dirac mass 
term ([7]) is even under time reversal T, the Haldane mass term ([5} is T-odd (for a recent review of the transformation 
properties of different mass terms in graphene, see Ref. [I(|). It is noticeable that the mass A was first discussed 
long ago in connection with inducing the Chern-Simons term in the effective action of (2 + l)-dimcnsional relativistic 
gauge field theories [3l| . 

III. QUASIPARTICLE PROPAGATOR 

The inverse bare quasiparticle propagator in the mixed (w, r)-representation is given by 

l S-V;r,r') = [(w + m)t° - imB-fo 3 + v F (n ■ j)]6(r - r'). (9) 

(To simplify the representation of all formulas, we omit the spin index in this section.) Similarly, the general structure 
of the full fcrmion propagator for quasiparticles of a fixed spin has the following form: 

iG- 1 (u;T,i f ) = |7°o; + w F F + (7r-7) + E + }(5(r-r'), (10) 

where F + and E + can be viewed as generalized wave-function renormalization and self-energy operators, respectively. 
In the special case of the bare propagator in Eq. ([5]), the corresponding functions are F^ are = 1 and S^ arc = (ji — 
(x B Ba 3 )^ . 

By definition, F + and S + arc functions of the three mutually commuting dimensionless operators: (tv ■ 7) 2 ^ 2 , 7° 
and is±7 1 7 2 , where sj_ = sgn(eB) and £ = \Jfic/\eB\ is the magnetic length. Taking into account that (7 ) 2 = 1 and 
(is^-y 1 ^ 2 ) 2 = 1, the operators F + and S + can be written in the following form: 

F+ = f + 7 °.g + Mj/yVff + ^7°7V/, (H) 
£+ = A + 7 // + is ± 7 1 7 2 M + is_ L 7°7 1 7 2 A, (12) 

where /, /, g, g, A, A, /i, and ft are functions of only one operator, (tt ■ ~/) 2 £ 2 . Note that for the functions /1, ft and 
A, A we keep the same notations as for the parameters fx, ft and A, A in Eqs. ([5]) and Eqs. ([Tj) , (jSJ), respectively. 
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It is obvious from the representations in Eqs. ([TO and (IT21) that F + and S + do not necessarily commute with 
(tt ■ 7). It is convenient, therefore, to introduce two other functions F~ and E~ , which satisfy the relations: 

^+(7^7) = (TVf)F-, (13) 
S+(tt-7) = (7r-7)S". (14) 

As follows from their definition, the explicit representations of these functions read: 

F- = / - 7 °<7 - is ±1 V.9 + is ±7 °7 V/, (15) 
XT = A — 7 // — is J _7 1 7 2 /i + is J _7°7 1 7 2 A. (16) 

These are obtained from F + and S + by reversing the signs in front of the two terms that anticommutc with (tv ■ 7). 

The physical meaning of the functions A, A, /i, and jl that appear in the definition of 5^ is straightforward: A 
is the Dirac mass function, A is the Haldane (time-reversal odd) mass function, fi is the charge density chemical 
potential, and ji is the chemical potential for charge density imbalance between the two valleys in the Brillouin zone. 
As for the functions /, /, g, and g that appear in the definition of F^ , they are various structures in the wave function 
renormalization operator. 

It may appear that, in the most general case, the full propagator (|10p can also include another wave-function 
renormalization, multiplying the frequency term 7°w. This is not the case, however, because this Dirac structure 
is already included in the self-energy E + , which may depend on uj in general. We note at the same time that the 
solution for E + will turn out to be independent of ui in the instantaneous approximation utilized in this study. This 
fact is also one of the reasons that makes it particularly convenient to separate the term 7°w from the generalized 
self-energy operator S + in (fT0|) . 

As mentioned earlier, functions /, /, g, g, A, A, /1, and jl are functions of (tv ■ ~/) 2 £ 2 , whose eigenvalues are 
nonpositive even integers: — 2n = —(2N + 1 + .s^si2), where N = 0, 1,2, . . . is the orbital quantum number and 
s\2 = ±1 is the sign of the pseudospin projection. Therefore, in what follows, it will be convenient to use the following 
eigenvalues of the operators F^ and (sec Appendix lAl for more details), 

F,^ ee fn + so9„ + s 12 g n + s s 12 f n , (17) 
S so, S i2 = + S /i„ + Si2/U„ + s s 12 A n , (18) 

where /„, /„, g n , g n , A„, A„, /i„, and jl n are the eigenvalues of the corresponding coefficient operators in the nth LL 
state. Further, so = ±1 and s\2 = ±1 are the eigenvalues of 7 and is±7 1 7 2 , respectively. 
In terms of eigenvalues, the inverse propagator is derived in Appendix [A} Its final form reads 

*G- 1 (w;r,r / ) = e^'^iG-^wjr-r 7 ), (19) 

iG-\lo-v) = ^E E E {«owin(f)+ [soMn.a + A^] [6%L n ($ + 5%L n ^(0] 



n—0 a— ±1 so — il 



+ ^(7 • r)(/„,a - *05n,«r)^-l(0}^o,.o«r. (20) 



where L" are Lagucrre polynomials (L„ = L n ). We also introduced the following short-hand notations: 

i=ijL 2p-' Hr,r') = - S J X + X '^~ y '\ (Schwinger phase) (21) 

and 

Mn,<r = A*n + vfin, A„ iCT = A„ + (7 A„ , (22) 

f 7l,(T = fn + f/n, 9n,a = 9n + <^9n- (23) 

Note that, by definition, the Laguerre polynomials L" with negative n are identically zero. Finally, V So , Sl2 are the 
projectors in the Dirac space, 

V So , Sl2 = ^(l + .so7o)(l + si2is_L7 1 7 2 ), with s ,si2=±l- (24) 
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Similarly, the expression for the propagator itself reads 

G(w;r/) = e 4 * (ry) G(w;r-r'), (25) 

— £/2 00 f 1 \ X 



ri=0<r=±l Sn=±l 



-pin • + » )2_ £a / n-i(0}^, (26) 



(w + m„, ct ) 2 - £ 2 „ 
where the energies in the lowest and higher LLs are 

Eo, a = ctA , ct = Ao + ctAq, (27) 



£ n , CT = y j2n(v 2 F /P)[fl„-gl <T }+Al rr , for n > 1. (28) 
The corresponding energies of quasiparticles are determined by the location of the poles of propagator (|26|) . i.e., 

w , CT = -Mo.o- + Eo,(T, (29) 

Un,o = -Mn, CT ± -E'n.er, for U > 1. (30) 

Let us note that c = ±1 is the eigenvalue of matrix isj_7 7 7 = 5^7375, which up to the overall sign s± is the 
quantum number associated with the valley. This follows from the explicit representation in Eq. ([6]) and from our 
convention for the four-component Dirac spinor, whose first two components are associated with valley K and the 
last two components with valley K' . 

IV. SCHWINGER-DYSON EQUATION 

The Schwinger-Dyson (gap) equation for the fermion propagator in the random-phase approximation (RPA) is 
shown diagrammatically in Fig. [TJ Note that in contrast to the naive mean-field approximation, the RPA Coulomb 
interaction includes the polarization (screening) effects, which are not negligible in the dynamics responsible for 
symmetry breaking in graphene. 

It is important to emphasize that the gap equation for the fermion propagator in Fig. [T] contains two tadpole 
diagrams. One of them is connected with the Hartree contribution due to dynamical charge carriers, while the other 
takes into account the background charge from the ions in graphene and in the substrate. The presence of both 
tadpoles is essential to insure the overall neutrality of the sample. Indeed, the equation for the gauge field implies 
that the two tadpole contributions should exactly cancel that yields (the Gauss law): 

j o ° xt - e Tr[ 7 G] =0. (31) 

As clear from the above arguments, this is directly related to the gauge symmetry in the model. (Since the Gauss 
law does not take place in models with contact interactions (27| . there is an analogue of only one tadpole diagram 
describing the Hartree interaction, which contributes to the gap equation.) Thus, the resulting Schwinger-Dyson 
equation for the fermion propagator takes the form 

G-\t - t'; r, r') = S'^t - t'; r, r') + e 2 7 ° G(t - t'; r, r') 7 °L>(i' - t; r' - r), (32) 

where G(t;r,r') is the full fermion propagator and D(t;r) is the propagator mediating the Coulomb interaction. 

At this point it is instructive to compare the cases of a local four-fermion interaction and a nonlocal Coulomb 
interaction. In the case of a local four-fermion interaction, the right-hand side of the Schwinger-Dyson equation 
contains <5(r — r') and the fermion propagator only at the point of coincidence G(0;r,r). This means that the right- 
hand side of the Schwinger-Dyson equation is a constant in the momentum space and, hence, it does not renormalize 
the kinetic tv ■ 7 part of the fermion propagator (i.e., F + = 1). Also, in the case of a local four fermion interaction, 
S + does not depend on the LL index n. Clearly, this simplifies a lot the analysis of the gap equation. The situation 
changes in the case of the nonlocal Coulomb interaction. 

Following Ref. [loj . we consider the instantanenous approximation for the Coulomb interaction by neglecting the 
dependence of the photon polarization function H(cj, k) on u. Then, in momentum space, the photon propagator 
takes the following form: 

D ^ k ^ D ^ = eo[ k + mk)y (33) 
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FIG. 1: The diagrammatic form of the Schwinger-Dyson equations for the electron and photon propagators in the mean field 
approximation. 



where 11(0, k) is the static polarization function and eo is a dielectric constant. In essence, the instantaneous approx- 
imation neglects the retardation of the interaction. This may be a reasonable approximation for graphene, whose 
charge carriers propagate much slower than the speed of light. It should be kept in mind, however, that such an 
approximation has a tendency to underestimate the strength of the Coulomb interaction (lfj . I32I . [33j . 

Unfortunately, it is difficult to find exact solutions of Eq. (|32[) . Therefore, one has to use some approximations. 
Here, we will study the dynamical symmetry breaking in the model under consideration retaining contributions only 
of the lowest and several first LLs. Obviously, this approximation is consistent only if the dynamically generated 
gaps are suppressed compared to the Landau scale Ei = y/hv F \eB\/c (it characterizes the energy spectrum of the free 
theory and is the gap between the lowest and first LLs). 

In the instantaneous approximation, the photon propagator reads 

d 2 k f duj n/ t^-iut+ik-r i f°° dk kJ (kr) 



By noting that the Schwinger phase on both sides of the gap equation (J32J) is the same, we arrive at the following gap 
equation for the translation invariant part of the propagator: 

■n-U \ -c-i/- ^ e ' r dn rdkkMkr) ~ 

iG (uj;r)=iS w;r / 75-/ r ^n/n u T G ii - r T ■ 35 

e J_ lx> 2iiJ 2?r k + 11(0, k) 

As shown in Appendix [B] this is equivalent to an infinite set of algebraic equations, see Eqs. (|B23I) . (|B24I) . (IB26[) and 
(TB271) . 



V. NUMERICAL RESULTS 



In this section, we present our numerical solutions to the truncated set of gap equations. In contrast to our previous 
analysis in Ref. [27j . where a model with a contact interaction was studied in detail, here we investigate the effect 
of the long-range interaction on the dynamics of symmetry breaking in QH effect in graphene. For this purpose, we 
will neglect the screening effects, captured by the polarization function 11(0, k). This approximation may be justified 
in the case of a strong magnetic field, which is of our main interest. In the case of fi = 0, for example, one has 
11(0, k) ~ ak(k£) in infrared [l(| and, therefore, the effects of the polarization function are suppressed with respect to 
the bare Coulomb interaction indeed. This point greatly simplifies the numerical analysis because now the kernel of 
the gap equation takes an analytical form. As for a more rigorous study of the dynamics with screening (e.g., in the 
random phase approximation), it will be presented elsewhere [34j|- Note that the strength of the Coulomb interaction 
is characterized by the graphenc's "fine structure constant" a = e 2 /c^vf, which is approximately equal to 2.2/eo- In 
our numerical calculations below, we use eo = 1. 

In order to use numerical calculations efficiently, it is important to understand all energy scales in the problem at 
hand. Ignoring the large energy cutoff due to a finite width of the conductance band, there are essentially only two 
characteristic energy scales in the action of graphene: (i) the Landau energy scale eg = y/hvp\eB\/c and (ii) the much 
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smaller Zccman energy Z = hbB. In our numerical calculations, we measure all physical quantities with the units of 
energy in units of eg. Numerically, these are 

eg = yJhv 2 F \eB\/c = 26^ B[T] meV, Z = \ibB = 5.8 x 10~ 2 B[T] meV, (36) 

where B[T] is the value of the magnetic field measured in Teslas. The corresponding temperature scales are Sg/ks = 
SOOy^Bpl K and Z/k B = 0.67B[T] K. 

Note that the magnetic length £ = y/hc/\eB\ and the Landau energy scale eg are related through the Fermi velocity 
as follows: t = hvp/eg = 26 ran/ yj B[T] , where we used vf/c = 1/300. 



A. Renormalization of the Fermi velocity (weak field) 

Let us start from the simplest analysis, when the role of the dynamical mass parameters is negligible and there is 
no significant splitting of the LLs. This is presumably the case when the magnetic field is not so strong. Even in this 
case, however, there is a very interesting dynamics responsible for the renormalization of the Fermi velocity. This 
is also interesting from experimental point of view because the renormalized value of the Fermi velocity parameter, 
which is also a function of the Landau index n, affects the energies of optical transitions (3^413 • 

To this end, let us consider only the subset of the gap equations, which involve the wave-function renormalization 
/„, see Eq. (|B19[) . Even this subset contains an infinite number of gap equations for each choice of spin, i.e., 

[1 - n F (E n > — n) — tif (E n ' + fj)] , for n > 1, (37) 



2 4^, nVS 



where we took into account the definition of the coefficients K n \ n in Eq. (|Bllj) and used the approximate expression for 

the LL energies: E n . a rj E n = ^/2negf n ,a, which are independent of the valley quantum number a = sosi2- Note that 
the spin index is also omitted, which is justified especially at weak fields. Here we use the notation np(x) = l/(e x + 1) 
for the Fermi-Dirac distribution function. 

Before analyzing the set of equation for the wave-function renormalization /„ numerically, it is instructive to note 
that the expression on the right hand side of Eq. (|3T|) is formally divergent. Indeed, by taking into account the 
asymptotes of the kernel coefficients k^'-i n-i as n ' ~ > °°> scc Eq. (|B13|) . we find that the sum over n' on the right 
hand side of Eq. (|37|) is logarithmically divergent. From quantum field theoretical point of view, of course, this is just 



an indication that the coupling constant a is also subject to a renormalization [381 ] . For our purposes in this study, 
however, we may simply assume that the sum over the Landau levels is finite. Indeed, in contrast to actual relativistic 
models, the effective action for quasiparticles of graphene is valid only at sufficiently low energies. Moreover, it is also 
clear that the energy width of the conducting band of graphene is finite. In fact, it can be shown that the formal 
value of the cutoff in the summation over the LL index n is approximately given by n max ~ 10 4 /B[T] (39|, where 
B[T] is the value of the magnetic field in Teslas. In our numerical calculations, wc will use a much smaller cutoff 
n max . For all practical purposes, when dealing with the observables in lowest few LLs, such a limitation has little 
effect on the qualitative and in most cases even quantitative results. Thus, in the rest of this subsection, we choose 
the value of the cutoff to be n max = 100. (We checked that the numerical results for /„ with the cutoffs n max = 50 
and n max = 150 are qualitatively the same. The values of /„ have a tendency to grow with increasing n max . It is 
understood, of course, that such a growth should be compensated by the renormalization of the coupling constant in 
a more refined approximation.) 

The effective Fermi velocity in the nth LL is determined by the following relation: %>F, n = fn^F, where the 
numerical values of the wave function renormalization are shown in Fig. [2] There we show many sets of the results 
which correspond to different values of the chemical potentials. The points are the actual data, while the lines 
connecting the points are shown for eye guiding the data for fixed values of the chemical potentials. The data on the 
top line corresponds to small values of the chemical potential, |^| < V^eg. The other lines correspond to the chemical 
potentials in the energy gaps between nth and (n + l)th LL (with n = 0, 1, 2, . . . from top to bottom). A part of the 
same data is also given in Table |U 

In general, the renormalized Fermi velocity vf.u is about 10% to 20% larger than its nonrenormalized value vf- 
These results seem to be somewhat smaller than the predictions in Ref. [4pj . One should keep in mind, however, that 
there are considerable uncertainties in the theoretical predictions for the renormalized values of the Fermi velocity. In 
part, these are associated with a relatively large value of the coupling constant in graphene and with the logarithmic 
running of the wave function renormalization itself. 
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TABLE I: Values of the wave function renormalization /„ for several values of the chemical potentials. 
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\fj,\ < \[2e t 


1.270 


1.243 


1.227 


1.214 


1.205 


1.197 


1.190 


V2ei < M < V±£e 


1.194 


1.224 


1.217 


1.208 


1.201 


1.194 


1.188 


\fiet < n < VEei 


1.193 


1.224 


1.217 


1.208 


1.201 


1.194 


1.188 


V&ee < fi < VEee 


1.166 


1.177 


1.200 


1.199 


1.194 


1.189 


1.184 


V8e e < n < VWe e 


1.150 


1.156 


1.165 


1.184 


1.185 


1.182 


1.179 


VWei < ^ < Vl2e e 


1.149 


1.155 


1.164 


1.182 


1.184 


1.182 


1.179 


VV2e e < fi < -JTIei 


1.138 


1.142 


1.148 


1.156 


1.172 


1.174 


1.173 



A convenient quantitative measure of the many-particle effects in the transition energies is given by the prefactors 
Cn,n> , which are introduced as deviations from the non-interacting carriers in graphene, 

AE n , n , = E n , ± E n = (V2ri ± V^M) e e + aeiC n ^ , (38) 

where, once again, a small Zeeman splitting of the energy levels is ignored. Note that both terms in (|3"8")l scale as 
~ \f~B~ and experimental data on infrared spectroscopy of LLs of graphene clearly confirm this behavior [35l - [37j . On 
the other hand, the dependence of the coefficients C„ i)m on the LL pair allows one to get information on many-body 
effects. By making use of our notation for the wave- function renormalization, we obtain 

■ 1), (39) 

when n ^ and n' ^ 0. (For transitions from the n = level and for transitions to the n! = level, the LLL never 
gives any contribution to the corresponding prefactors.) By making use of our results for /„, we obtain the values of 
prefactors C n .n' ■ For transitions between several low-lying LLs, the values of the prefactors are listed in Table ILT1 



TABLE II: Values of some prefactors C n y , assuming that [i — 0. 



C-i,o 


C-ia 


C-2,0 


C_2,l 


C_2,2 


0.176 


0.352 


0.224 


0.400 


0.448 



It is clear that the Coulomb interaction contribution to the LL transitions slightly increases the transition energies 
above their noninteracting values in accordance with experimental data (e.g., see Fig. 3a in Ref. [36j]). 
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B. Lowest Landau level approximation (strong field) 



In the case of a strong magnetic field, the interaction of charge carriers plays a very important role. One of the 
empirical features in this regime is a lifted degeneracy of the Landau sublcvels, associated with spin and valley quantum 
numbers. This observation is naturally interpreted as an indication of symmetry breaking in the corresponding ground 
states with hxed filling factors. Here we analyze such a dynamics, using a low-energy effective model with a long-range 
Coulomb interaction in detail. 

To start with the analysis, let us first consider the solutions to the gap equations for a fixed spin of charge carrier. 
Keeping only the LLL terms (n = 0) in the infinite series on the right hand side of Eqs. (|B23[) and (|B24[) . we obtain 



M = 



"0,0 



Xeff 
^0 



ae e (o) 
o K o,o 



1 - n F (Kf - ) - n F (Af 



Mo 



(40) 
(41) 



where we used the notation ^Q ff = /j,q — Aq and Aq S = Aq — fto for the two independent combination of parameters 
that determine the spectrum of quasiparticle at the LLL sublevels (see Eq. (|29[1 ). 



XeS 

^0 ! 



-Mf 



Af. 



(42) 



At zero temperature, the gap equations become 



Mo 



IA 



off 



aee 

4v / 2tt' 



gn^oXK 



eff 



OLEi 
?=' 

4v / 2tt 



, eff 



1^1), 



IA 



off 



(43) 
(44) 



where we took into account the exact value of Kqq = l/(2-\/27r), see Table [TTTJ One of the solutions to this set of 
equations is of the "magnetic catalysis" type, i.e., 



Mo 



|Af | = 



/i, for 

as £ 
4V2T 



(45) 
(46) 



(Note that, according to this simplified estimate in the LLL approximation, the value of the gap |AQ ff | ranges from 
30 K to 190 K when the magnetic field changes from 1 T to 50 T, assuming that the coupling constant a is of order 
1.) The other two solutions are 



-co < n < 



4V2tt 

< [L < OO, 



Mf 



Mo 



M 

/' - 



ae t 



4V2tt' 
aee 



4V2tt 



A 



eff 



A 



off 



0. 
0. 



(47) 
(48) 



All of these solutions can be easily found numerically. In addition, the numerical study of the finite temperature 
equations shows two extra (probably unstable) solutions, satisfying the approximate condition ^g ff sa |Ag ff |. These 
solutions seem to survive even when one approaches the limit T — > 0. They are lost in the analytical study of the above 
zero-temperature equations because they correspond to the vanishing value of the argument in the step-functions. 
The ambiguity of the step-function in this case prevents us from finding the same solutions. Note that all 5 solutions 
can be continuously continued into each other through a sort of hysteresis loop without any discontinuities in the 
values of yff and |Ag ff |. This fact alone strongly indicates that some "intermediate" solutions are mctastable or even 
unstable. 



C. Numerical solutions for /i near the lowest Landau level 



Let us now perform a more realistic numerical analysis by including several LLs and accounting for the spin degree 
of freedom in the analysis. The explicit form of the gap equations is presented in Appendix [B] see Eqs. (|B23[) through 
(|B27[) . Because of the large number of the dynamical parameter, in this section, we will use a rather small number 
for the cutoff n max = 5 in the summation over the LLs. 
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-0.76 



-0.78 



-0.05 



0.00 
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FIG. 3: Numerical results for free energies of several lowest energy solutions, realized when the Fermi energy is near n = LL. 
The dashed lines correspond to hybrid solutions. 



The gap equations for both spins look the same except for the value of the chemical potential: it is ^ = fi — Z 
and fii = fi + Z. By repeating the same analysis as above, we find that there exist many more solutions around the 
vanishing value of fi. Keeping only a subset of several qualitatively different solutions with lowest energies, we find 
among them a pure Dirac mass solution, four types of the Haldane mass solutions and four types of hybrid solutions 
(see below). The pure Dirac mass solution has nonzero Dirac masses for both spins and no Haldane masses [i.e., the 
order parameters are triplets with respect to both SU^(2) and SU±(2)]. Four Haldane mass solutions are determined 
by four possible sign combinations of the two time-reversal breaking masses: (i) Ao,-|- > and A04 > 0; (ii) Ao,^ > 
and A04 < 0; (iii) A 0) -f < and A04 > 0; (iv) A ,^ < and A ,j. < 0. All of these are characterized by singlet order 
parameters with respect to both SU^(2) and SUi(2) symmetry groups. Similarly, four hybrid solutions are determined 
by the following conditions (i) A 0j -|- ^ and A 0i j, > 0; (ii) A ^ ^ and A 0,4, < 0; (iii) A ^ > and A 0i 4_ 7^ 0; (iv) 
Ao.t < and A04 ^ 0. The common feature of the hybrid solutions is that one of their order parameters is a triplet 
with respect to SU^{2) or SUi(2) group, while the other order parameter is a singlet with respect to the other group. 

Including many "intermediate" branches of solutions, in our numerical analysis we were able to identify several 
dozens of non-equivalent solutions in the vicinity of the LLL. The free energies for all of these solutions can be easily 
calculated using Eq. (|C5|) in Appendix [C] 

The free energies of several lowest energy solutions are plotted in Fig. [3l Four singlet type solutions and one triplet 
solution are shown by solid lines in the figure. Four lowest energy hybrid solution are shown by dashed lines. There 
are three qualitatively different regions, in which the lowest energy states are different, i.e., 



At the points fi = ±Z , there also exist hybrid solutions with the same lowest values of the energy as the two Haldane 
mass solutions. 

• v = —2 (the LLL is empty). This is a singlet solution, whose free energy as a function of [i is shown by the 
green solid line in Fig.[3j This solution corresponds to the unbroken U^(2) x J7j.(2) symmetry and has the lowest 
energy for n < —Z: 



M < -Z 
Z <n<Z 

^ > z 



A , T > & A 04 > 0, 
A 0iT > & A04 < 0, 
A , T < & A04 < 0. 



(49) 
(50) 
(51) 




0, M g ff t = M - Z - Ao, tl A , t w 0.225e,, 
0, n* s i = n + Z-A , i , A 04 « 0.225£,. 



(52) 
(53) 



In this case the LLL quasiparticle energies are 



t*>o,t = + Z + |Ao,f| > 0, 
wo,4. = -M--Z+IA04I > 0, 



(x2), 
(x2). 



(54) 
(55) 



Since none of the LLL sublcvcls arc occupied, this solution corresponds to a v = 



—2 state. 
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v — — 1 (the LLL is one-quater filled). This is a hybrid solution, whose free energy as a function of \i is shown 
by the red dashed line in Fig. [3J The symmetry in the corresponding state is spontaneously broken down to 

Uf(2) x (1) x U^ K where the two latter factors describe U(l) transformations at a fixed spin and a 
fixed valley. At \x = — Z, this solution is degenerate in energy with the solutions for the v = — 2 and v = cases 
and is given by 

Ao ff T = 0, ^ ff t = n - Z - A 0lt , Ao.fw 0.225e* (56) 
AgJJ, « 0.225e*, ^ = /i + Z, A 04 = 0. (57) 

In this case the LLL quasiparticle energies arc 

w ,t = + |A , t | > 0, (x2), (58) 
wo,; = -M--^ + Ao4>0, (59) 
W04 = -n-Z-h^<Q. (60) 
Since only one LLL sublevel is occupied, this solution corresponds to a v = — 1 state. 

v = (the half filled LLL; neutral point). This is a singlet solution. Its free energy as a function of fj, is shown 
by the blue solid line in Fig. [3] The symmetry in the corresponding state is Uf(2) x f/j.(2), but the Zeeman 
splitting is dynamically enhanced. This solution has the lowest energy for — Z < fx < Z 

AEf t = 0, Mat = I 1 - Z - A °,t> A o,T « 0.225e £ , (61) 
A<u = 0, /iS4 = M + ^-Ao,t, A 04 « -0.225^. (62) 
In this case the LLL quasiparticle energies are 

w ,t = + |A , t | > 0, (x2), (63) 
w<U = -M-^-|A ^| <0, (x2). (64) 
Since two of the LLL sublevels are occupied, this solution corresponds to a v = state. 

v — 1 (the LLL is three-quater filled). Similarly to the solution for the v — —\ state, this is a hybrid solution. 
Its free energy as a function of /j, is shown by the gray dashed line in Fig. [3] The symmetry in the corresponding 
state is spontaneously broken down to lli K \l) x lli K \l) x U±{2). At fx = Z, this solution is degenerate in 
energy with the solutions for the v = and v = 2 cases and is given by 

Ajft « 0.225e £ , ^ = n - Z, A i)t = 0, (65) 

A^ = 0, ^ ff t = M + ^-A ,t, A 04 w -0.225e<. (66) 

In this case the LLL quasiparticle energies are 

wo,f = -A* + 2 + Ajf t > 0, (67) 

w , t = -n + Z- A£ ff t < 0, (68) 

wcu = -M-^-|AoJ <0, (x2). (69) 
Since three of the LLL sublevels are occupied, this solution corresponds to a v = 1 state. 

z/ = 2 (the LLL is filled). This is another singlet solution. Its free energy as a function of /i is shown by the 

light-brown solid line in Fig. [3l This solution corresponds to the unbroken f7f(2) x U\.{2) symmetry and has the 
lowest energy for \i > Z: 

A° ff t = 0, ^ = fi - Z - A ,t, A , t ~ -0.225e £ , (70) 

A ct = 0, Lil s i = Li + Z-A , i , A 04 « -0.225e*. (71) 

In this case the LLL quasiparticle energies are 

w 0lt = -M + Z-|A 0jt | <0, (x2), (72) 
wo.4. = -M-^-|Ao4|<0, (x2). (73) 

Since all LLL sublevels are occupied, this solution corresponds to a v = 2 state. 
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D. Numerical solutions in the n = 1 Landau level 

In this subsection, we present the lowest energy solutions in the n = 1 LL, when this level is partially or completely 
filled (see Fig. @|. In this case, we were also able to identify several dozen non-equivalent branches of solutions. 
Using Eq. (|C5[) . we selected the ones with the lowest free energies. The main features of the corresponding solutions 
are described below. Note that filling of the n = 1 LL leads to changing the dispersion relations also in other LLs. 
Therefore, although only parameters A 1; A 1; 0* and fi determine the QH plateaus at the n = 1 LL, we write down 
their values for two neighbor levels, the LLL and the n = 2 LL. (Recall that, according to Eqs. ([T7| and (|18|) . A„, 
A„, /i„, /„ are functions of the LL index n.) Physically, the information concerning the gaps in other LLs is relevant 
for experiments connected with transitions between Landau levels [36| . 

• v — 3 (the n — 1 LL is one-quater filled). Hybrid solution, which is valid for ~ V2~££ — Z . The symmetry in 



this state is spontaneously broken down to Uf(2) x u[ K \l) x u[ K 

Af f = 0, Vo S t = Mt ^ A o,t, A 0lt = -0.225s* (74) 

A^ = 0.052s* l i% s i = » i -A , i , A 04 = -0.277s* (75) 

A 1)t = 0, 1)t = A*t + 0.053s* A 1)t = -0.067s* / 1)t = 1.142, (76) 

A u = -0.018s* 014 = 04 + 0.148s* A 14 = -0.049s* / u = 1.103, (77) 

A a ,t = 0, M2,t = A*t + 0.040^, A 2 , t = -0.051s* / 2 , t = 1.111, (78) 

A 24 = -0.006s* A*24 = H + 0.091s* A 24 = -0.045s* f M = 1.102. (79) 

• 1/ = 4 (the n = 1 LL is half filled). Singlet solution, which is valid for \[2ei — Z < \i < V2sf + While 
formally the symmetry of this state is the same as in the action, Uf(2) x U±(2), it is characterized by a dynamically 
enhanced Zeeman splitting: 

AS ff t = 0, ^ ff t = ^t-Ao4, A , t = -0.225s* (80) 

A°o4-0, M c 04 = 04 - A04, Ao4 = -0.328s* (81) 

A14 = 0, 0i, t = /if + 0.053s* A Xit = -0.067s* f 1A = 1.142, (82) 

Ai4 = 0, 0i4 = 04 + 0.244s* A14 = -0.031s* /14 = 1.065 (83) 

A24 = 0, 024 = 0t + 0.040s* A24 = -0.051s* / 2 4 = 1.111, (84) 

A 2 4 = 0, 024 = 04 + 0.142e* A 24 = -0.039s* f %i = 1.092. (85) 

• v — 5(the n = 1 LL is three-quater filled). Hybrid solution, which is valid for ~ \/2ei + Z . The symmetry in 
this state is spontaneously broken down to U^ K \l) x U^ K \l) x Ui(2): 

Ag ff f = 0.052s* 0ff T = t -Ao4, Aq4 = -0.277s* (86) 

AS ff 4 = 0, 0g ff 4 = 04-Ao4, Aq4 = -0.328s* (87) 

A14 = -0.018s* 014 = t + 0.148s* A14 = -0.049s* /14 = 1.103, (88) 

Ai4 = 0, 014 = 04 + 0.244s £ , A 14 = -0.031s* /14 = 1.065, (89) 

A24 = -0.006s* 024 = 0t + 0.091s* A24 = -0.045s* / 2)t = 1.102, (90) 

A 2 4 = 0, 024 = 04 + 0.142s* A 2 4 = -0.039s* ^4 = 1.092. (91) 

• v = 6 (the n = 1 LL is filled). Singlet solution with the unbroken U^(2) x £/|,(2) symmetry. It is valid for 
> v2s^ + Z.that was mentioned in the previous subsection. 

Ajf t = 0, 0jf t = t - A04, Aq4 = -0.328s* (92) 

A°o ff 4 = 0, 0° o ff 4 = 04 - A04, Ao4 = -0.328s* (93) 

A14 = 0, 014 = t + 0.244s* A14 = -0.031s* /14 = 1.065, (94) 

Ai4 = 0, U = 04+ 0.244s* A u = -0.031s* /14 = 1.065, (95) 

A 2 4 = 0, 024 = t + 0.142s* A24 = -0.039s* f 2tt = 1.092, (96) 

A 2 4 = 0, 024 = 04 + 0.142s* A 2 4 = -0.039s* / 24 = 1.092. (97) 
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1.30 1.35 1.40 1.45 1.50 1.55 

FIG. 4: Numerical results for the free energies of several lowest energy solutions, realized when the Fermi energy is near n = 1 
LL. 



Before concluding this section, it is appropriate to mention that, while all the results for the QH states, associated 
with filling the n = 1 LLs, arc qualitatively similar to those obtained in Ref. [27]], the Coulomb long-range interaction 
makes all gaps and other dynamical parameters functions of the LL index n. 



VI. DISCUSSION 



The present analysis of integer QH plateaus connected with lifting the fourfold degeneracy of LLs leads to the 
following results. The v = QH state is a QH ferromagnet with significantly dynamically enhanced Zeeman splitting. 
The v = 1 QH plateau is described by a CDW for the electrons with a spin polarized along the magnetic field. Such a 
CDW is absent for the electrons with the opposite polarization. The solution corresponding to the v = 2 QH plateau 
is a singlet and it necessarily includes the Haldane mass. At the n = 1 LL, the dynamically generated gaps are much 
less than those at the LLL. Therefore, the corresponding QH states could be observed only in rather strong magnetic 
fields and very clean samples. 

It is noteworthy that the present results obtained in the mean-field approximation with the long-range Coulomb 
interaction are qualitatively similar to those obtained in Ref. 27j , where a short-range four-fermion interaction was 
used. The main reasons of this fact are: (a) the universality of the magnetic catalysis phenomenon, and (b) in both 
these models, there is essentially one dimensionfull parameter, the Landau scale eg = yjTiv 2 F \eB\/ c, through which the 
energy spectra of the QH dynamics are expressed. In connection with that, we recall that the dimensionfull coupling 
constant Gi nt in the model in Ref. [13] is taken as G mt ~ A/e^, and the free dimensionless coupling A plays the role of 
the Coulomb coupling there. On the other hand, it is noticeable that while the gaps arc constant for the short-range 
interaction, they decrease with increasing the LL index n in the case of the long-range Coulomb one. Note that 
without magnetic field energy gaps are functions of a frequency u) and a momentum k for long-range interactions fiol ] . 
The discrete index n replaces the momentum when magnetic field is switched on. 

It should be noted that a recent study of the cyclotron mass in suspended graphene [4l[ suggests a very strong 
renormalization of the Fermi velocity (at the electron density n = 10 9 cm -2 it is almost three times larger than 
the value commonly used). We note however that the effect of the Fermi velocity renormalization was studied in 
several other experiments, such as photoemission spectroscopy [42| . scanning single-electron transistor measurements 
of the local electronic compressibility 43 1, and optical infrared measurements in a magnetic field [3514371 [44[. While 
the optical measurements observed a departure from the dispersion for nonintcracting particles, its value was small. 
The pre sent study also suggests that the corresponding renormalization effects are much smaller than those found in 
Ref. [4lj . and, in fact, they are comparable to those observed in Refs. [36l. Eij]. 

It is instructive to compare the present results with the realization of the magnetic catalysis in (3 + l)-dimcnsional 
relativistic theories, such as QED [45j] and QCD (46[. A dynamical generation of only one type of the Dirac mass 
(m = A) is usually considered in relativistic theories in 3 + 1 dimensions. As we know, the set of the order parameters 
in graphene is much richer. It is only recently that new order parameters, analogous to the QHF order parameters and 
the Haldane mass, have been studied in QED and QCD [47], HH]. In large part, the studies of the QHE in graphene 
have been responsible for those new studies in relativistic theories. [Note, however, that in 3 + 1 dimensions the analog 
of the Haldane mass term describes an axial- vector current density, rather than a mass [48j .] 
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Appendix A: Fermion Green's function 



1. General structure 



The general structure of the inverse full fermion propagator (for a single spin species) in a magnetic field is written 
in Eq. (fT0|). Making use of the two sets of operators, F^ 1 and X^, defined in Eqs. (fTTj) . (fl"2j) . (fi"5|) , (fl"6j) . we derive the 
following formal representation for the full Green's function: 

-l 



G = 



j°uj + v f F+(-k -7) + S H 



j°uj + v f F+(-k -7) -XT | 7°w + vfF + (tt • 7) + S h 
~/°lj + v f F + (-k -7) -XT |cj 2 + 7°cj (jh + - XT 



7°w + v f F + (tt -7) - XT 

„ „ A . "I -1 

F+F-(ir-j) 2 -E+XT j 



f 



where 



7 °u;(X] + -XT) = 
F + F~ = 

E+XT = 



2uj((i + isj_7°7 1 7 2 /i), 
(/ + mx/tV/) 2 - (.9 + ^i 7 °7 1 7 2 g) 2 
f + ~f -9 2 -~9 2 + 2*.s ± 7°7 1 7 2 (// - 99), 
(A + is ± 7°7 1 7 2 A) 2 -(fi + isj.7 7 1 7 2 A) 2 
A 2 + A 2 - y 2 - fi 2 + 2i.s ± 7°7 1 7 2 (AA - nji), 



(tt • 7) 2 r = -7r 2 r - (e/c)Bi7 1 7 = -ttV - is_L7 7 , 



(Al) 

(A2) 

(A3) 

(A4) 
(A5) 



and we recall the definition of the magnetic length i = ^jTic/\eB\ and sj_ = sgn(eB). We also used the fact that the 
operator 7T 2 ^ 2 has the eigenvalues (2N + 1) with N = 0, 1, 2, ... . In the Landau gauge, A = (0, Bx), utilized here, 
the corresponding normalized wave functions read 



where Hj^(x) are the Hermite polynomials. These wave functions satisfy the conditions of normalizability 

d 2 rtl}* Np {r)ijj Nlp >(Y) = 5 NN 'S(p-p'), 



and completeness 



00 

00 

/ dpiP Np (r)r Np (r') = 5(r-r'). 

N=0 „ 



(A6) 



(A7) 



(A8) 



Note that operators F^ and S 1 * 1 , when acting on eigenstates \N,p, so, S12), reduce to their eigenvalues: 



F 



-(2iV + 1 + s 12 ), s , flu] \N,p, s , S12) 



F± [(7r-7) 2 £ 2 ,7o,«s±7V] |iV,J5, s , S12) 

= ^+(' S t 2 S +l)/2l Ar 'P' S 0' S 12>, 

S± [(7r-7) 2 £ 2 , 7o ,is ± 7 1 7 2 ] \N,p,s ,s 12 ) = ±± [-{2N + 1 + s 12 ), s , s 12 \ \N,p,s ,s 12 ) 

= ^'£^ 1)/9 \N iP ,s 0i a 12 ). 



(A9) 
(A10) 
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Taking into account that si 2 = ±1, we conclude that the LL index n = N + (,si2 + l)/2 is an nonnegative integer. 
Therefore, as follows from the definitions in Eqs. pT|) . (|T2"j). ([15)) and (jTHJ) , the eigenvalues are determined by the 
following expressions: 



F°°' S12 = f n + s g n + s 12 g n + s s 12 f n , 
E^ 12 ee A„ + s /x„ + s 12 /t„ + s s 12 A r , 



(All) 
(A12) 



where /„,/„, 9ni 9n, A„, A n , /i„, and jl n are the eigenvalues of the corresponding coefficient operators in the nth LL 
state. 



2. Inverse propagator 

Let us introduce the following projectors in the Dirac space: 

V So ,s 12 = ^(1 + so7o)(l + si2«s J _7 1 7 2 ), with s , s 12 = ±1, 
Then, we can write the inverse fermion propagator as 

iG-V; r,r') = £ i G-^ ia (w; r, r')7>. ,. ia) 

S ,Si 2 =±l 

(note that the projector is multiplied on the right and that the ordering is important) where 

00 „ 

e -£/2+»*(r,r') 



n=0 



2t^ 2 


z 






g-£/2+iS(r 








g-£/2+iS(r 
j : 





aow + S^" ( 1 / la+1)/2 



MO 

£(l + Sl2 >;WCrV(r-rO^(0 
n=0 
00 

X) (1 - si2)^^;°i +1 7 • (r - r')M(0 



n'=0 



+ 



e -€/2+i$(r,r') 1 + 
-f/2+i*(r,r') 

2tt£ 2 



00 ,• 

X + Kti 1 } M0 + ^7 ■ (r - rOiCr^O 



n=0 



X j [«o* + E*"" 1 ] MO + ^ 7 ■ (r - r')F- + s ?' +1 in(0 



n=0 



(A13) 
(A14) 



-£/2+i*(r.r') 00 r 

= ^ E + 5- 12 iK°> Sl2 + S^KiT] MO + -f 7 ■ (r - t>)K$- m L\S& \ , (A15) 

with 72/ = n — 1 (the prime is omitted after the first use) and the short-hand notation 

(r-rT , {x + x'){y-y') 

? 2 ^2 . n r , r J 2^ 2 

In the derivation, we changed the integration variable p — > q = p£ + — j Sj _£zJL anc j took into account that 



(A16) 



e -^r(.x+p£ 2 ) 2 e -^z(x'+pe 2 ) 2 e ™±p(y-y') = e -g 2 -£/2-M*(r,r') 



(7T • 7)-0ATp 
7T 2 ljj Np 



I 

2f 

'11 

2( 
i 

2? 



y/2(N + TjipN+i, P - V2N^ N - hp 
y/2(N + l)i>N+l,p + V2Nif> N - llP 



(A17) 
(A18) 
(A19) 



y / 2(N + l)(j 1 - is ±l2 )4> N+1 . p - V2N(ji + is±72)^JV-i lP ] ,(A20) 
2N + 1 



(A21) 
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and integrated over the quantum number p by making use of the formula 7.378 in Ref. |50| . 

r x2 H m (x + y)H n (x + z)dx = 2 ,l 7r 1 / 2 m!z"- m J L^- m (-2yz), (A22) 



assuming m < n. Here are the generalized Laguerre polynomials, and L n = L? n . 

The final representation of the inverse propagator is in the form of the Schwinger phase and the translation invariant 
part, 

iG-\u;;r,r') = e^'^i G~V; r - r'), (A23) 

— f /2 00 



n— s ,s 12 — ±1 



-^(7 • r)F-;'r^Ll(0}Vs ,s 12 , (A24) 



or, it can be written in alternative form as Eq. (|20j) in the main text. The inverse of the bare propagator has a similar 
structure: 

iS' 1 (uj;r,r') = e^AirV;!--!-'), (A25) 

— £/2 00 

iS-\w,T) = E E {«OwMO + *OM[*i^«(0+^»-x(0] 



n— (T— ±1 so— ±1 



£(7 ■*)£»-! (Oj^W, (A26) 



It should be emphasized that the translation invariant part G (ui; r) is not the inverse of G(w; r) used later. 



3. Propagator 



Using the identity in Eq. (|A1[) . the spectral expansion of the unit operator (|A8[) and the projectors in Eq. (|A13I) . 
we can write the propagator as follows: 



G(w;r,r') = ^ G So , Sl2 (w; r, r')P So , Sl2 , 

S ,Sl2=±l 



(A27) 



where 



+ (1 + S12)/2,S0S12 



,-f/2+t*(r,r') 



2nt 2 

,-£/2+»*(r,r') 



OO 

E [ S ° W ~ S n+(Ir7 a +\)/2 



n=0 

00 

E(l-^Xf7-(r-r 



A^ n+ (l +Sl2 )/2,s s 12 



n'=0 



-e/2+i*(r,r') j + 
27Tf 2 2~ 



00 



v — so, — 1 

sow - E n+1 ' 



111 

r- 



Mn+l, ao 



e -e/2+i*(r,r') 1 



S12 



2tt^ 2 



M n ,-s 



■ (r - r ')F~ S0 ' +1 



(A28) 
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where 

M n ,a = {uJ + H n ,a? with <? = ±1, (A29) 

E , a = a Ao,cr = Aq + a Ao, (A30) 



V = VM^)[/^-S^]+^. for »>!■ (A31) 
The energy spectrum at LLL and higher LL's is determined by the location of the poles in Eq. (|A28|) , 

Wo, CT = -Mo,<r + Eo,a, (A32) 
ut,a = _ Mn,<r ± ^W.cr, for 71 > 1. (A33) 

The final representation of the propagator can be written in the form of the Schwinger phase and the translation 
invariant part, 

G(u;r,r') = e l<E> ( r < r '>G(u;; r - r'), (A34) 

— £/2 OO / \ A 

<W> - *^E E E { y + ^r^" [^(0 + ^Ln-M 

+ ^-'' (^^> - K) K- (A35) 

Note that the free propagator reads 
S(w;r,r') = e l ^ ry *> S{tu; r - r'), (A36) 

Ofu-ri : - e ' €/2 V V r «o(" + /4[*iWO+*ffWC)] fag, ^-!(Q 1 

^ ' J ~ 2t^ 2 ^ ^ I (w + m) 2 - 2m; 2 /£ 2 + p W ' r ^ (w + y _ 2nv z ,p J 7 «o,«» 

2tt£ 2 ^l 7 (w + /i) 2 -2m; 2 /£ 2 £2 a ; ^ + ^ 2 _ 2m; 2 /i 2 /' { ' 

where s\2 = ctsq and V± = (1 ± is± r y 1 j 2 )/2. 



Appendix B: Gap equations 



The general form of the gap equation is given in Eq. (|32p . By making use of the photon propagator (|34[) . describing 
the Coulomb interaction in the instantaneous approximation, and the mixed (u), r)-representations of the fermion 
propagators, we arrive at the following gap equation: 



*<?">; r)=»S->;r)-- / — 



dn r°° dk kj (kr) ~ 



2tt J 2nk + U(0,k) 



7 u G(C!;r) 7 u 



(Bl) 



Multiplying both sides of the gap equation (|B1[) by either e £/ 2 £ n (£) or e ^ 2 (7 • r)i^(^), and then integrating over 
r, we find that this is equivalent to following set of equations: 



Mn,cr — M — oA„ i0 



S s ° + 



„2 00 r°° 



iwv-/* + ^^ = ~L/ 3 ^y 2 ^ + n(o,fc) 



dfi f°° dk kC { °) n (kl) 



f2 + fj, 



2 oo „oo 



n — 



i fA n / CT - f2 + /Z n '+1 a- — cA n /_)_i CT 
_„ H 



fn,a + Sa9n,o ~ 1 



dfi /"» dfc *4^-l,n-l( W ) /«> + SOS* 



neo £± 1 J_ 00 2n Jo 2tt fc + n(0,fc) M„ v 



(B2) 

for n > 1 , (B3) 



where 



1 r 2 



2£ 2 



^" [ ^ ) ^ ( TO ) W - (-l) m+ "e-^ 



fc 2 ^ 5 



fc 2 ^ 2 



(B4) 
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= 



2? 4 Jo 



drr e ^ L„ 



1 ^) L 1 

2i 2 J n \2P 



Jo(fcr) = (-ir+"(m + 1) e"— 

= (-ir + "(n + l)e-^Lr 



k 2 £ 2 
k 2 l 2 



t m — n 
L n 

j m—t 



,2p2- 



2 



.(B5) 



To obtain the results on the right hand sides, we used the table integral 7.422 2 in Ref. [5C 
The gap equations can be equivalently rewritten as follows: 



Mn,<7 - M - a ^n,a 
Pn,a — M + oh n%a 



E(0) f°° dfl 
K n'n / 2= 



dflil + p n <. a + crA„/ CT 



n'=0 

oo 



E( ) Z" 00 dft 

K n'-l,n-l / 2 7T 
i'=l J -°° Z7F 



■M„', CT 
<ir2 + jJL n i a — CrAn/^ 



for n > 0, 



for ?i > 1, 



/n,<r = 1 + ta£f 2^ _ / 



_ x 2tt M n < 



, for n > 1, 



,9n.cr = «ae<? 



E — 



where a = e 2 /(eo«F), = ^f/^> and 



,(p) 



, 2tt -M„' )(T ' 
m. n (hi) 



for n. > 1, 



P = 0,1. 



(B6) 
(B7) 
(B8) 
(B9) 

(BIO) 



l 2vr k + 11(0, A;) 

When the screening effects are neglected, i.e., n(0, k) = 0, we can use the explicit form for of Cm] n (with p = 0, 1) in 
Eqs. (|B4|) and (|B5[) and obtain the following analytical expressions for Km] n (with p = 0, 1): 



m.n 



oo oo 

dx 



n^o J 2ir 

T(p + l/2)r(l/2 + n)r(l/2 + to) 



2tt 
o o 



(-1) 



m+n 



2V2 



2 v / 27T 2 TO!n! 

nin(m,n) 

E 

fc=0 



3 F 2 (-to, -n, p + 1/2; 1/2 - to, 1/2 - n; 1) 
IV+l/2 + fc) 



(m - fc)!(n - fc)!r(l/2 - to + fc)r(l/2 - n + k)k\ 



(Bll) 



where we used formula 2.19.14.15 from Ref. [5l|. The values of Km] n (with p = 0, 1) at small values of to and n are 
given in Tables IIIII and IIVI The leading asymptotes for n — > 00 (at finite to) are 



,(o) 



2to - 1 



O 



n^o ~ 2Try/2n 87r(2n) 3 / 2 \n 5 / 2 

(m + 1) (m+l)(3m-l) 
n^o ~~ 4?rV2n 167r(2n) 3 / 2 



for n — > 00, 
for 



,5/2 



(B12) 
(B13) 



The zero temperature gap equations (|B6|) through (|B9|) arc straightforwardly generalized to the case of non-zero 
temperature by making the replacement il — > iQ m = i-nT{2m + 1) and using the Matsubara sums instead of the 
frequency integrations, 



dn 



£<...> ^ E <■■■)■ 



m— — 00 



Then, we use the following table sums: 

00 

t E 



* E 



1 1 — np(a + p) — np(a — p) 

(O m — i/u) 2 + a 2 2a 

iQ m + p up (a + p) — np(a — p) 



(Cl m -ip) 2 + a 2 



(B14) 

(B15) 
(B16) 
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TABLE TTL 

±iLiJi_JJ_J JL JL JL * 


v oiucii ui <"Tn.n 


when the effects of pofarization tensor are 


nptrlppf prl 




ri, m,n 


m = 


m = 1 


m = 2 


m = 3 


m = 4 


m = 5 


n = 
n = 1 
n = 2 
71 = 3 

rc. = 4 
n = 5 






3 


5 








4V2vr 


16V27T 

7 


32V2ir 
11 


256v / 2tt 


512V2tt 


4%/27T 

3 


H v"2- 
7 


32V2ir 
41 


64V2ir 
51 


512%/27r 
329 


1024\/27r 
569 


16V2tt 
5 


32V2tt 
11 


128V2ir 
51 


256V2tt 
147 


2048V2ir 
759 


4096V27T 
1245 


32V2ir 
35 


64V2n- 
75 


256V2ir 
329 


512V2n- 
759 


4096V2jt 
8649 


8192V2vr 
11445 


256^2tt 
63 


512\/2jr 
133 


2048v / 2tt 
569 


4096^2tt 
1245 


32768^2jt 
11445 


65536\/27r 
32307 


512V2vr 


1024V2ir 


4096\/27r 


8192V27T 


65536-/2-7T 


131072 \Z2tt' 



TABLE IV: Values of Km,n when the effects of polarization tensor are neglected 





771 = 


m = 1 


m = 2 


771 = 3 


771 = 4 


771 = 5 


n = 


i 


i 


3 


5 


35 


63 


4v/2tt 


8^27T 


32V2if 


64^211 


512^211 


1024 v / 2tt 


n = 1 


1 


7 


15 


23 


155 


273 


&V2w 


16s/2ir 


64%/2tt 


128V2ir 


1024V2ir 


2048V2tt 


n = 2 


3 


15 


153 


171 


1065 


1809 


32^2if 


64V2~7r 


256V2vr 


512^2jr 


4096^2if 


8192^2tt 


n = 3 


5 


23 


171 


759 


3495 


5505 


64V2ir 


128-/27T 


512V2vr 


1024\/27r 


8192V2ir 


16384\/27r 


n = 4 


35 


155 


1065 


3495 


57225 


67365 


512^211 


1024^271 


4096 V2tt 


8192^2tt 


65536^2tt 


131072^211 


72 = 5 


63 


273 


1809 


5505 


67365 


261207 


1024V27T 


2048V2tt 


8192V2n- 


16384V2tt 


131072 v^x 


262144 v^-ir 



and derive the finite temperature gap equations, 



(l>n,a — A* — cA n a = 



fJ«n,,c? — M + crA n , c 



oo 

^ ^n'.n^F {En', a - AV.cr) - »*F {E n > ,cr + Mn',<r) 
ra'=0 

" ,(T [1 - HF (-En'.a - AV,<t) ~ »F (-E'rt'.a + \, for 71 > 0, 

tin', a > 
oo 

X! K n°'-X,n-l{ n -F (^',(7 ~ A*n',a) - ™F (-En'.a + AV.ct) 

1 - n F (-B„', CT - M„',a) - n F (E n ^ a + n n ', a )) }, for n > 1, 



(B17) 



+ 



crA, 



/n,cr 1 



oo (1) , 



(B18) 



E 



n E r , 



[1 - n_F (i?„' i(T - /i„', CT ) - n F {E n i >(T + Hn',a)\ , forn > 1,(B19) 



9n.<7 



oo (1) 

aee - — - A - 



2 f( n S„ 

n' — 1 



3 "^fL [1 _ nF (_E n , >ff _ ^ n , i<7 ) - n F {E n , t „ + Hn>,a)] , for n > 1. (B20) 



Note that the v! = term in the first equation can be rewritten in a simpler form: 



OlEl (0) 

2 



_ (0) , 

= ^ _K 0,7l tan 



2T 



2n F f a-Ao )(T - Mo,<rJ - 1 
(As defined in Eq. (|A30|) . E , a = fAo, CT Separating the LLL from the higher LLs, we derive 

oo 



(B2i; 



2T 



E n i. 



[1 - ?l F (J5 n / )<T - /!„>) - Tljr + /Un'.tr)] |, 



(B22) 
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which is equivalent to the following set of gap equations for the LLL parameters: 



-J- ^ K "°!o I"? C^n'^- _ Mn',+) - "-F + /Lt n ',+) + n F (-£„',_ - - ™F (-EW',- + AV,- 



E 



aee ^-^ 

n'=l 



A, 



4 ^ 



K n',0 



E n ' 



[1 - n F (£„',+ - Mn',+) - «F (-£«',+ + + 



[1 - n F (E n i- - /X„',_) - TIF {E n ',- + Mn',-)] i 



(B23) 



, eff 



[l - n F (Af - M S ff ) - n F (Af + 

CO 

~~T E ^,0 t nF C^™',- - Mn',-) - + Mn',-) ^ "-F ~ AV,+) + «F (-Ew,+ + Mn', + )] 



53 ^OTT^ t 1 - n F [E n ',+ - Mn',+) - "-F (£«',+ + Mn', + ) 
4 T-T. 



~F 53 K «"!o-^7- I 1 " "F (-En',- - Mn' -) - ™F (E n > _ + Mn' -)] 



(B24) 



where the two independent combinations of the LLL parameters in the gap equations are 

uf =fx - A , Af = A - Ac 



(B25) 



It is appropriate to mention that the parameters an d A cannot be unambiguously defined, while their combination 
/if can be and it is the only combination with a well defined physical meaning. The same is true for A and jj,Q, 
which are related to physical observables only through the effective Dirac mass Ag ff . 

The remaining equations for the chemical potentials and masses in the higher (n > 1) LLs read 



_ OLEl (o) 

Un,a — M + ^ K 0,n 



2n F (aAf- fJl f)-l 



oo (0) (0) 

Oi£t sr-^ K n'.n "+" K n'-l,n-l r , „ v , „ >, 



n'=l 



ae f « n / 



.(0) (0) 



^E 



n ^n' — l,n— 1 ^n'er 



A„, CT = [l - 2n F (<rAf - uf 



aet 



E 



„(0) (0) 
'V'.n ™n'-l,n-l 



«F {E n i^ a — Hn',a) — n F (E n ',a + Mn' 



(B26) 



00 ,c( ) 4_ J ) A 



n' = l 



En 



(B27) 



Appendix C: Free Energy Density 



Here we calculate the free energy density f2 in monolayer graphene with dynamically generated self-energy cor- 
rections and wave function renormalization by following the same approach, based on the Baym-Kadanoff formalism 
I52I1. that was used in Ref. [23 in the case of a model with a contact four-fcrmion interaction. In the model with the 
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r(G,D) = 






FIG. 5: The diagrammatic form for the effective action at two loops. 



Coulomb interaction, studied in this paper, the corresponding effective action T takes the form: 



T(G,D) = -iTr[LnG- 1 + S- 1 G-l]+-[LnD- 1 +D^ ) 1 D-l] 
- l JL [ d 3 u j d 3 u'tT[j°G(u,u')^°G(u',u)]D(u-u'), 



(CI) 



where u — (t,r). The trace, the logarithm, and the product S~ 1 G arc taken in the functional sense, and G = 
diag(G|, Gi). The diagrammatic form of this effective action is shown in Fig. [5] It is instructive to compare this 
action with the analogous action in the model with a contact four-fermion interaction [27J ■ In contrast to that model, 
here we have only one type of a two-loop diagram. It is responsible for the exchange interaction. The other "dumbbell" 
diagram (i.e., two fermion loops connected by the free photon propagator), which would be responsible for the direct 
(Hartree) interaction, does not appear in the present model with the gauge interaction. The absence of such a diagram 
is a formal consequence of the Gauss's neutrality condition, 



j c ° x tH-etr[ 7 G(u, W )] =0 , 



(C2) 



where Jc X t( w ) is the charge density of the substrate and carbon ions that compensate the overall nonzero charge 
density of free carriers in the graphene monolayer. 

The free energy density il is expressed through T as follows: SI = — T/TV, where TV is a space-time volume. When 
the fermion propagator satisfies the gap equation (|52")h the expression for the effective action T simplifies, i.e., 



-i Tr 



LnG- 1 + - (S^G-l) 



- [LmD- 1 + D-'D - 1] . 



(C3) 



Up to the one-loop photon contribution, which we will neglect in the following, this is formally the same result as in 
Ref. (27j . Therefore, we can use a similar derivation and change only the form of the Green's function G s (u,u') in 
order to take into account that a wave function renormalization is nontrivial here and that all dynamical parameters 
are functions of the Landau index. We finally obtain 



n = - 



1 



4tt£ 2 

OO 



^ \ Sgn(/X ,er) (M0,<r + - CrA 



10, cr 



<7=±1 k 



lMo,o 



|A 0>CT | 
2n/ n , 



En,, 



Sgn(A 0itT ) A 0:Cr - cr (p + (J, 0j<r ) 0(|Ao,er| - \(Mj,tr\) 



(C4) 



\E nta \ - \Hn,a\) - 2V2ne e 



where we introduced the function 9(n — 1), which is defined so that 6{n — 1) = for n = and 8(n— 1) = 1 for n > 1. 
At a non-zero temperature, we make the replacement u> — > iuj m = i?rT(2m + 1) and use Matsubara sums instead of 
the frequency integrations, and obtain 







4 ^ Yl \ + » - a ^°>") tanh 



a=±l 



M0,cr — crAo,c 

2T 



! ^2 [(Mn,<7 + /- 1 ) [n>F{En,o - Mn,<r) ~ np{E n ,cr + A'ri.a)] - 2 \/2nE e 1 - 2n F (V2ne e ) 
2n/„ iff e 2 



E n,a 



E n , 



[1 - n F (E rha - fin >a .) - n F (E n ^ + \i riy 



(C5) 
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